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ABSTRACT 

We have performed a high mass and force resolution simulation of an idealized galaxy forming from 
dissipational collapse of gas embedded in a spherical dark matter halo. The simulation includes star 
formation and effects of stellar feedback. In our simulation a stellar disk forms with a surface density 
profile consisting of an inner exponential breaking to a steeper outer exponential. The break forms 
early on and persists throughout the evolution, moving outward as more gas is able to cool and add 
mass to the disk. The parameters of the break are in excellent agreement with observations. The 
break corresponds with a rapid drop in the star formation rate associated with a drop in the cooled 
gas surface density, but the outer exponential is populated by stars that were scattered outward on 
nearly circular orbits from the inner disk by spiral arms. The resulting profile and its associated break 
are therefore a consequence of the interplay between a radial star formation cutoff and redistribution 
of stellar mass by secular processes. A consequence of such evolution is a sharp change in the radial 
mean stellar age profile at the break radius. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: structure — galaxies: pho- 
tometry — galaxies: spiral — stellar dynamics 



1. INTRODUCTION 



Since the early work of Ide Vaucouleurd (|1958l ) it 
has been recognized that the disks of spiral galax- 
ies generally follow an exponential radial surface 
brightness profile, and various theories have explored 
the physi cal causes and consequences of this prop- 
erty (e.g. iFall fc Efstathioul [19801 : iLin fc Pringlel fl987l: 
Dalcanton et alj Il997t iMo et all 11998b Ivan den Boschl 
since Ivan der KruitJ (|1979l Il987fi ~ 
that the outer regions of disks ex- 



disks provide a direct view of disk assembly. 

Several the ories for the formatio n of breaks have been 
investigated. IVan der KruitJ (|1987[ ) proposed that angu- 
lar momentum conservation in a collapsing, uniformly 
rotating cloud naturally gives rise to disk brea ks at 
roughly 4.5 scale radii. IVan den Boschl ([2001) sug- 



2001). However, 
has been known 



it 



hibit more varied behavior. This ha s been confirmed by 
an abundance of recent data (e.g. , Pohlen et al. 120001 
2002t iBland-Hawthorn et alj l2005t lErwin et all 120051 : 
Pohlen fc Truiilloll2006l hereafter PT06). In a sample of 



gested that breaks are due to angular momentum 
cutoffs of the cooled gas. More commonly breaks 
have been attributed to a threshold for star forma- 
tion (SF), wheth er due to l ow ga s density which sta- 
bilizes the disk (jKennicuttl Il989f). or to a lack of a 



cool equilibrium IS M phase ( Elmegreen fc Parravanol 
1994 ISchavd 12004). Us ing a semi-analytic model, 



Elmegreen fc Hunterl (|2006f) demonstrate that a double 



nearby late- type galaxies from the Sloan Digital Sky Sur- 
vey, PT06 found that about 60% have an inner exponen- 
tial followed by a steeper outer exponential (downward- 
bending), ^30% have the inner exponential followed by 
a shallower outer exponential (upward-bending), while 
only ~ 10% have no detectable breaks. Therefore breaks 
are a common feature of disk galaxies that any com- 
plete theory of galaxy formation must be able to explain. 
Furthermore, the discovery of U V emission at radii well 
beyond the Ha emission cutoff (|Gil de Paz et all 120051 : 
iThilker et al.ll2005l 120071). the observational eviden ce for 
inside-out disk growth (jMunoz-Mateos et al.ll2007f), and 
detec t ions of disk breaks in th e distant universe (jPerezI 
12004 iTruiillo fc Pohlenl I2005D . suggest that the outer 
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exponential profile may result from a multi-component 
star formation prescription. The e xistence of extended 
UV disks (e.g IThilker et all I2007D and the lack of a 
clear correlation of Ha cut-offs and optical dis k breaks 
(|Pohlen et all 12004 iHunter fc Elmegreen! l2006f) further 
complicates the picture. Regardless, while a sharp SF 
cutoff may explain the disk truncation, it does not pro- 
vide a compelling explanation for ex tended outer ex- 
ponen tial components. Alternatively, iDebattista et all 
(2006) demonstrated that the redistribution of angular 
momentum by spirals during bar formation also produces 
realistic breaks in collisionless V-body simulations. 

In this Letter we present the first results from a se- 
ries of high-resolution Smooth Particle Hydrodynamics 
(SPH) simulations of isolated galaxy formation aimed 
at exploring the formation and evolution of disk breaks 
and outskirts in a massive, high surface brightness galaxy 
without a strong central bar. Resulting breaks are anal- 
ogous to downward-bending breaks seen in observations. 
The clear advantage of our approach over past attempts 
is that we use a fully self-consistent physical model of the 
system, making no a priori assumptions about the distri- 
bution of material in the disk. Rather, we allow the disk 
to grow spontaneously under the effects of gravity and 
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gas hydrodynamics, itself influenced by star formation 
and feedback. The iV-body approach (at sufficiently high 
mass and force resolution) ensures that we capture the 
dynamical processes contributing to disk evolution. Fur- 
thermore, the inclusion of prescriptions for SF and feed- 
back allows us to make observational predictions across 
the break region. We show that (1) the break forms 
rapidly (< lGyr) and persists throughout the evolution 
of the system, moving outward as the disk mass grows; 
(2) the break is seeded by a sharp decrease in star forma- 
tion which is caused in our simulation by a rapid decrease 
in the surface density of cool gas; (3) the outer disk is 
populated by stars that have migrated, on nearly circular 
orbits, from the inner disk, and consequently the break 
is associated with a sharp change in the radial mean stel- 
lar age profile; (4) break parameters agree with current 
observations. 

2. SIMULATION METHODOLOGY 

Resolved stellar population data in disk outskirts, 
which are now becoming available (e.g. Ferguson et al.l 
[20061: Barker et all [20071 : Ide Jong et all 120071 ). provide 
strong constraints on theories of break formation. There- 
fore, the inclusion of SF and feedback is required to as- 
sess break formation models. For this reason we have 
run simulations of gas cooling, collapsing, and forming 
stars inside a live dark matter halo within which the gas 
is initially in hydrostatic equilibrium. This has the fur- 
ther advantage of making no assumptions about the an- 
gular momentum distribution within th e disk, which can 
stron gly affect its subsequent evolution (Ibebattist a et al.l 
[20061) . 

We construct initial conditions as in [Kaufmann et all 
(|2007f h The initial system c onsists of a virialized spheri- 
cal NFW dark matter halo ([Navarro et al.|[l997l ) and an 
embedded spherical hot baryonic component containing 
10% of the total mass and following the same density dis- 
tribution, which at the end of the simulation yields a disk 
mass fraction of ~ 5%. The mass within the virial radius 
is 10 12 M©. A temperature gradient in the gas compo- 
nent ensures an initial gas pressure equilibrium for an 
adiabatic equation of state. Velocities of gas particles 
are initialized according to a cosmologically-motivated 
specific angular momentum distribution with jar and 
an overall spin parameter A = (j /G)^\E\/M 3 = 0.039 
(|Bullock et al.|[200ll ). Each component is modeled by 10 6 
particles; the dark matter halo is composed of two shells, 
with the inner halo of 9 x 10 5 particles (each of mass 
~ 1 x 10 6 Mq) extending to 200 kpc and an outer halo 
of 1 x 10 5 particles (each of mass ~ 3.5 x 10 6 M Q ) be- 
yond. All gas particles initially have a mass of 1.4 x 10 5 
Mq. We use a softening length of 100 pc for all dark 
matter particles and 50 pc for baryonic par t icles. We 
adopt the best-fit values from IStinson et alj (|2006l ) for 
the parameters governing the physics of star formation 
and feedback. Our cooling prescriptions do not account 
for effects of UV background or metal line cooling. The 
global criteria for SF are that a gas particle has to have 
n > 0.1 cm~ 3 , T < 15000 K and be a part of a converging 
flow; efficiency of star formation is 0.05, i.e. 5% of gas 
eligible to form stars is converted into stars per dynami- 
cal time. Star particles form with an initial mass of 1/3 
gas particle mass, which at our resolution corresponds to 
4.6 x 10 4 Mq. A gas particle may spawn multiple star 




Fig. 1. — Azimuthaly-averaged properties of the disk 1, 5, and 
10 Gyr. The top panels show the stellar surface density profiles, 
with the dotted lines representing double exponential fits. The 
two exponential components were fit simultaneously to the profile, 
excluding the innermost and outermost regions. The point of in- 
tersection of the two exponentials is taken as the break radius. In 
all panels, the vertical lines indicate the location of the break. The 
second row shows the surface density of cool gas. In the third row, 
we show the SFR density calculated from the stellar mass formed 
in the previous 10 Myr. The mean stellar age (normalized by time 
of output) as a function of radius is shown in the bottom row. 

particles but to avoid gas particles of unreasonably small 
mass, the minimum gas particle mass is restricted to 1/5 
of the initial mass. The simul ation is evolved with the 
parallel SPH code gasoline (jWadslev et al.l 12001 ) for 
10 Gyr. 

We have also performed simulations with 10 5 and 10 7 
particles per component, thus bracketing our fiducial 
run. While the details of the gas cooling are somewhat 
resolution-dependent ([Kaufmann et al.l l2006| ) resulting 
in morphological differences, we find convergence in the 
modeling of overall halo cooling, star formation and its 
dependence on gas surface density, and disk dynamics. 
We therefore consider the conclusions and predictions of 
this Letter numerically robust. The 10 6 particle reso- 
lution used here represents a compromise between ade- 
quate statistics for detailed analysis in the outer disk and 
computational cost for 10 Gyr of evolution. 

Our simulation should be thought of as modeling disk 
formation and evolution from the cooling of hot gas after 
the last gas-rich m ajor merger, as suggested by cosmo- 
logical simulations ([Brook et al.ll2004 120071 ). Although 
we make use of simplifications such as an initially spher- 
ical distribution of matter, lack of halo substructure, 
the ignoring of subsequent accretion of dark matter and 
baryons, and an initial gas density profile that mimics 
that of the dark matter, studying the idealized isolated 
case allows us to analyze in detail the important dynami- 
cal processes affecting the evolution of a massive isolated 
disk. The lessons we learn from this idealized case will 
later be applied to galaxies evolved in a full cosmological 
context (R. Roskar et al. in preparation). 

3. BREAK FORMATION AND THE OUTER DISK 

In Figure [TJ the break is already evident as soon as a 
stable disk forms at 1 Gyr, moving outward as the disk 
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Fig. 2. — 2D histogram of final particle radii vs. particle for- 
mation radii for particles on nearly circular orbits at 10 Gyr. The 
right panel shows only particles beyond the break. 

grows and persisting throughout the simulation. A sharp 
drop in the local SFR is always present at the break ra- 
dius. The drop in SFR is not due to a volume density 
threshold for star formation, but is instead associated 
with a rapid decrease in the gas surface density (the star 
formation follows a Kennicutt-Schmidt law at all times) 
and a corresponding sharp increase in the Toomre Q pa- 
rameter. We verified that the break is not seeded by our 
star formation recipe by running several simulations with 
different values of the threshold density and found that 
the location of the break did not depend on the particu- 
lar value used. Since density is inversely proportional to 
radius, the cooling time increases outward. By construc- 
tion, the angular momentum is directly proportional to 
cylindrical radius, which means that higher angular mo- 
mentum material will take longer to cool. This leads to 
the radial extent of the gas disk (and subsequently the 
star-forming disk) being limited by the maximum angu- 
lar momentum of material that has been able to cool. 

Figure [1] shows clearly that a sharp change in the mean 
stellar age profile coincides with the location of the break. 
We now investigate the cause of this ubiquitous feature. 
Figure [2] shows a 2D mass-weighted histogram of final 
(t = 1Q Gyr) versus formation radii of stellar particles on 
nearly circular orbits, defined as having J Z /J C (E) > 0.9, 
where J C (E) is the angular momentum of a circular or- 
bit with the energy of a given particle. This criterion is 
satisfied by ~ 80% of the outer disk particles. Therefore 
the outer disk is a rotationally-supported structure. It 
is clear from Figure [2] that significant redistribution has 
occurred throughout the disk; including all of the parti- 
cles in this figure, the root-mean-square change in radius 
during the lifetime of each particle [(Ai?) 2 ] 1 / 2 = 2.4 kpc 
(mean Ai? ~ 0). Most strikingly, nearly all (~ 85%) par- 
ticles beyond the break (Figure [21 right) have migrated 
there from the inner disk. The mean epicycle radius of 
these outer disk particles is only ~ 2 kpc, while the mean 
Ai? = 3.7 kpc. Hence, the presence in the outer disk 
of particles which formed in the inner disk cannot be 
attributed to radial heating alone but is a signature of 
s ecular radial redistri b ution . 

Scll wood fc Binnevl (|2002l ) describe how transient spi- 
ral arms move corotating stars radially without heating 
the disk. Fourier analysis shows that spirals of a wide 
range of pattern speeds are present in different parts of 
the disk at all times, enabling shuffling of stars at all 
radii. As the disk size increases, spirals of ever larger 
corotation radii and lower pattern speeds become com- 
mon, enabling transport of material over large distances. 
The largest migrations coincide with the growth of strong 
spirals, confirming that the migrating stars are in reso- 



nance with the spirals. Thus radial migration is common 
and stars of all ages are continuously transported to well 
outside the break. Details of these mechanisms will be 
presented in a future paper. 

Since the surface density break also corresponds to a 
break in SFR, the outer disk is relatively deficient in 
young stars. The outermost regions of the disk are pop- 
ulated by old stars which migrated from the inner disk, 
since they have had the most time to reach there. This 
redistribution naturally leads to the minimum in stellar 
age at the break and a positive age gradient in the outer 
disk. Although the importance of spiral arms for the evo- 
lu tion of the outer disk h as already been demonstrated 
bv lDebattista et all (|2006f ). in their models the break was 
formed by resonant coupling between the central bar and 
spiral arms. The cause of the break in our model is fun- 
damentally different since it is ultimately caused by a 
drop in the SFR, while the outer exponential is popu- 
lated by secularly redistributed stars. While a recurring 
oval does a ppear in the disk, a stro ng bar like the ones 
described in lDebattista et aLl (|2006f ) does not form here. 

The radial migrations are not directed only outward; 
indeed spirals cause much more radial mixing than trans- 
port. FigureOshows that the radial changes are symmet- 
ric for particles with small deviations Ai? and at smaller 
radii. However, for particles with larger Ai? or larger 
initial radii, more particles move outward than inward. 
Part of the excess angular momentum required for this 
transport is provided by the gas, which is funneled to- 
ward the center. 

4. COMPARISON WITH OBSERVATIONS 

Our simulation is not intended to reproduce all prop- 
erties of real galaxies since we lack the full cosmological 
context. Nevertheless, the evolution of disks in isolated 
galaxies since the epoch of last major mergers at z ~ 2 
(i.e. about 10 Gyr) should not be very different from 
that in our model. We may therefore anticipate that our 
model matches nature well if it captures the essential 
physics. 

We have investigated how the break structural pa- 
rameters of our simulated disk at different points in 
the simulatio n comp are to those observed by PT06 and 
iPohlen et al.l (|2002r) . Ratios of break radius to inner 
scalelength agree well with the data from PT06 through- 
out the simulation. The mean ratio of break radius to 
inner scalelength, R{> r /hi ni in our simulation is 2.6, while 
PT06 find a mean value of 2.5, an indication that our sim- 
u lation captures the co rrect scale of the phenomenon. 

Ide Jong et al.l (|2007f) presented radial star count pro- 
files of the edge-on galaxy NGC 4244 for different stellar 
populations. Their main result was that the break occurs 
at the same radius regardless of age and distance from 
the mid-plane. In Figure [3] we show edge-on density pro- 
files for stars from our mode l with age bins rough ly cor- 
responding to those used by Ide Jong et al.l (|2007f) . Solid 
lines correspond to stars in the midplane (\z\ < 0.5 kpc), 
dashed lines to stars at 0.5 < \z\ < 1.0 kpc, and dotted 
lines to stars at 1.0 < \z\ < 1.5 kpc. The vertical line 
marks the break as determined from the face-on profile. 
The break occurs at the same radius for all age bins, but 
becomes less pronounced with age due to radial redistri- 
bution. In the two oldest age bins (where we have enough 
particles to make a meaningful comparison) , the break is 
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Fig. 3. — Stellar density profiles for stars in different age bins. 
An arbitrary offset has been applied for clarity. See text for details. 

also constant with height away from the midplane. Tran- 
sient spirals are responsible for the largest amounts of 
redistribution, which is very efficient and can operate on 
relatively short timescales. Therefore the agreement of 
the location of the break between stars of different ages 
is a natural consequence of migration caused by transient 
spirals in a growing disk. 

The decrease in mean stellar age with radius is unsur- 
prising if disks form inside-out, and such trends have in- 
deed been observed (Ide Jong 1996; Bell & dc Jong l2000l ; 
iMacArthur et all 12004 iMunoz-Mateos et all 120071 ). 
Thus, the negative age gradient out to the break is ex- 
pected, but the abrupt change in the gradient and its 
close correlation with the break radius has not previously 
been noted. However, a positive age gradient beyond the 
break has been observed in M33 (B arker et al.ll2007t ). al- 
though, to the best of our knowledge, no complementary 
age profile exists for the inner disk of that galaxy. 
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5. CONCLUSIONS 

We have shown that in a self-consistent model, where 
the stellar disk forms through gas cooling and subsequent 
star formation within a dark matter halo, breaks in the 
stellar surface density form through the combination of 
different effects. A rapid drop in the SFR seeds the break 
and secular evolution populates the outer exponential. In 
our model the SFR drop is due wholly to a drop in the 
surface density of gas. However a break in SFR induced 
by other means (e.g. a volume density threshold or per- 
haps warps) would lead to similar behavior of the outer 
disk and stellar density break parameters. Our model 
properties satisfy current observational constraints, both 
in the statistical sense of break properties from galaxies 
in SDSS, as well as the much more detailed observations 
of breaks in stellar populations of NGC 4244. Though 
our model does not account for the effects of evolution in 
a full cosmological context, its simplicity assures that this 
is the minimal degree of evolution (with no interactions 
or formation of a significant bar) and should therefore 
be rather generic provided that a disk is massive enough 
for strong transient spirals to form. The model predicts 
that there should be an abrupt change in the radial mean 
stellar age profile coincident with the break, which can 
be tested with future observations. 
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